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We present results of LDA calculations (band structure, densities of states, Fermi surfaces) for possible 
iron based superconductor BaFe2Se3 (Bal23) in normal (paramagnetic) phase. Results are briefly compared 
with similar data on prototype BaFe2As2 and (K,Cs)Fe2Se2 superconductors. Without doping this system is 
antiferromagnetic with ~250K and rather complicated magnetic structure. Neutron diffraction experi- 

ments indicated the possibility of two possible spin structures (antiferromagnetically ordered "plaquettes" or 
"zigzags"), indistinguishable by neutron scattering. Using LSDA calculated exchange parameters we estimate 
Neel temperatures for both spin structures within the molecular field approximation and show n ("plaque- 
ttes") spin configuration to be more favorable than T2 ("zigzags"). 

PACS: 71.20.-b, 74.70.-b, 75.10.-b, 75.25.+Z 
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1. INTRODUCTION 

During the recent years novel iron based high- 
temperature superconductors [1] attracted a huge num- 
ber of both experimental and theoretical investigations 
(for review see [2, 3]). This flow of papers is still continu- 
ing uninterrupted. Major chemical classes of iron based 
superconductors are pnictides [1] and chalcogenides [4] 
(for comparative study see Ref. [5]). 

Recently BaFe2Se3 (Bal23) system [6] was proposed 
'as possible superconductor (analogous to (K,Cs)Fe2Se2) 
iwith experimental indications of Tc ~11K. However, 
'later study claimed no superconductivity in this system 
down to 1.8K [7]. Both papers [6, 7] exploited neu- 
tron diffraction and discovered Bal23 to be antiferro- 
magnetic spin-ladder system with Neel temperature of 
about 250K, and proposed two possible spin configu- 
rations. In Ref. [6] these spin configurations (antifer- 
romagnetically ordered "plaquettes" or "zigzags" ) were 
reported to be equivalent in a sense of diffraction pic- 
ture and thus indistinguishable. Ref. [7] is rather in 
favour of "plaquettes" pattern. 

Most recent investigations of the Bal23 system 
Ishown that with doping (Fc deficiency) it becomes a 
semiconductor with the energy gap about 0.2 eV at 
room temperature [8, 9]. In Ref. [9] LSDA calculations 
were also used to calculate total energies of phases with 
different magnetic structures. From these results it can 
be seen that more favorable is apparently the "plaque- 
ttes" configuration. 
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In this work we present LDA calculated elec- 
tronic structure, densities of states and Fermi surfaces 
for Bal23 system in (high-temperature) paramagnetic 
phase (probably most promising for superconductivity) 
and compare the results with those for the previously 
studied systems, such as BaFe2As2 and (K,Cs)Fe2Se2. 
For antiferromagnetic phase we use LSDA calculated 
values of Heisenberg model exchange parameters to cal- 
culate Neel temperatures for "plaquettes" and "zigzags" 
spin structures, to compare relative stability of both 
structures. To this end we conclude that "plaquettes" 
configuration is more favorable. 

2. ELECTRONIC STRUCTURE 

We start from crystal structure of Bal23, which be- 
long to orthorhombic Pnma space group [6] . Main struc- 
tural motiff of the Bal23 system is a two-leg ladder 
going along 6- axis. This ladder is formed by Fe ions 
which are surrounded by tetrahedrally coordinated Se 
ions. Those ladders form a kind of "checker board" 
structure perpendicular to fe-axis. This structure is obvi- 
ously quite different from previously studied BaFe2As2 
and AFe2Se2 systems (see Refs. [5, 10, 11]), which have 
body centered tetragonal crystal structure. 

Using this experimentally established crystal struc- 
ture (lattice parameters and all atomic positions used in 
our calculations were taken from Ref. [6] ) we performed 
LDA band structure calculations within the linearized 
muffin-tin orbitals method (LMTO) [12] using default 
settings. 

In Fig. 1 we present LDA calculated band disper- 
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sions (on the right side) and densities of states (DOS) 
(on the left side). Similarly to Bal22 [10] and AFe2Se2 
[11] systems electronic states on the Fermi level and 
around it (from 1 to -2 eV) are mostly formed by Fe-3d 
orbitals. The Se-4p orbitals form bands at energies be- 
low -2 eV. Hybridization between Fe-3d and Se-4p states 
is rather moderate. 

DOS, states/eV/cell 
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Fig. 1. LDA calculated band dispersions and densities 
of states of paramagnetic Bal23. The Fermi level Ef 
is at zero energy. 

Further on in Fig. 2 we show band dispersions of 
Bal23 in close proximity of the Fermi level. We can 
see that these band dispersions are very much different 
from typical bands of previously studied pnictides or 
chalcogenides [10, 11]. Close to F-point there are two 
electron pockets and near the Brillouin zone border (Y- 
point), there are three hole-like pockets (see also Fig. 3, 
lower panel). Close to the Fermi level there are several 
band crossings and Van-Hove singularities. Thus Fermi 
surface topology can be changed rather easy upon dop- 
ing as in AFe2Se2 system [11], stressing the potential 
importance of doping for superconductivity search. 
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Fig. 2. LDA calculated band dispersions in the vicinity 
of the Fermi level for Bal23. The Fermi level is at zero 
energy. 



In Fig. 3 we present LDA calculated Fermi surfaces 
(FS) of Bal23. Overall shape of the Bal23 FS (upper 
panel of Fig. 4) is very different from those in typical 
iron pnictides or chalcogenides [10, 11]. First of all it 
looks pretty three-dimensional and does not have any 
well established cylinders. Near F-point we see a jabot- 
like structure and near Y-point we observe a kind of 
"bottomless pot" FS sheet. In the lower panel of Fig. 3 
we show FS cross-section for the — plane. Fermi 
sheets closest to F-point are electron like, while those 
FS sheets near Y-point are hole-like. 




Fig. 3. LDA calculated Fermi surface in the Brillouin 
zone (upper panel) and fcz = plane Fermi surface 
crossection (lower panel). 



3. MAGNETIC STRUCTURE 



Neutron scattering experiments [6, 7] on magnetic 
structure of Bal23 at T < T^^^^ ~250K can be inter- 
preted as indicating two possible types of spin config- 
urations, corresponding to irreducible representations 
Ti and T2 of the Pnma space group, giving the same 
diffraction peaks structure. Since neutron diffraction 
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Ti spin configuration of Bal23 ( "plaquettes" ) 
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T2 spin configuration of Bal23 ("zigzags") 
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Fig. 4. Experimental variants of AFM spin configura- 
tions Ti (upper panel) and T2 (lower panel) with LSDA 
calculated Hciscubcrg exchange parameter values (K). 
"Spin up" orientation corresponds to circles and "spin 
down" to squares. Figures in circles and squares enu- 
merate magnetic sublattuces. Dashed lines separate dif- 
ferent ladders. 



In case of magnetic structure with several coUinear 
magnetic sublattices Heisenberg model Hamiltonian for 
classical spin vectors can be written as: 



H 



'2 ^ ^ Ii,n,j,mSi,n^j,m ' 



(1) 



where Ii^n,j,m is exchange integral between site i of sub- 
lattice n and site j of sublattice m, S - classical spin 
vector length and ej,„ is the unit vector of spin direction. 
Self-consistent equations for thermodynamic averages of 
spin 2:-projections f7i„ =< > within mean-field ap- 
proach can be linearized near Neel temperature Tjv and 
written as: 



S 



(2) 



Due to translation invariance of the crystal the values 
of CTj„ = (Jn are independent of site number in magnetic 
sublattice. Then Neel temperature Tjv is determined by 
the solution of the full system of linear equations for (7„ 
for all magnetic sublattices. 

The values of exchange integrals for both spin config- 
urations (ti and T2) in Bal23 antiferromagnet for clas- 
sical Heisenberg model at T = were calculated using 
the method proposed in Ref. [13]. The results are shown 
by numbers (in K) on bonds in Fig. 4. 

Note that both n and T2 structures can be repre- 
sented as a set of 16 collinear sublattices with 8 with 
"spin up" and 8 with "spin down" orientations. The 
absolute value of 5 at T = on each Fe ion in all sub- 
lattices is the same. However spins on different sublat- 
tices arc connected with their magnetic surroundings by 
different sets of exchange integrals. Thus, at finite tem- 



perature the values |(T„ 



< Si„ > I will be different 



for different sites. We rcstri(;t ourselves to the account 
of only nearest and next nearest exchange bonds of a 
given ion in the spin-ladder where it belongs and also 
coupling of nearest neighbors in adjacent ladders. 

For Ti spin structure, using the symmetry of nearest 
magnetic environment (see Fig. 4, upper panel), we can 
see that the following relations should be fulfilled: 



experiments can not distinguish between these magnetic ci = CT3 = — (T5 = —ay = — aio = —<Ji2 = <yiA = fig 

structures, it is tempting to find more stable configu- (j2= 04, = —cfq — ~<Js = —erg = — cth = ais = (T15. (3) 
ration theoretically. One can do that by computation 

of Neel temperature for both configurations within the Then the linearized equation for the average value of z- 

molecular field approximation, using LSDA calculated projection of cti spin on the magnetic sublattice number 

Heisenberg exchange interaction parameters. 1 becomes: 
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+Jl,40-4 + Ji.sfg + 2Ji 14CT14 + 2Ji_i2<Ti2 + 2Ji_i6Cri6] ,(4) 

where Ji 2: and jj'y are nearest exchange parame- 
ters for selected site of the sublattice 1, with neighbors 
in 2, 3 and 7 sublattices within the ladder, to which 
the selected site belongs; Ji^ and Ji^g are next nearest 
exchange parameters within the ladder to which the se- 
lected site belongs; Ji 14, Ji,i2 and Ji,i6 are adjacent 
ladder nearest exchange parameters. 

Using Eq. (3) the Eq. (4) can be reduced to: 

TjvCTi = — [( Ji^3 + Jij + 2Ji 14 — 2Ji^i2 + 2Ji^i6)cri 

+ {Jl2 + J1A-Jl,s)<72]- (5) 

Similar equation for a2 takes the form: 

Tn<72 = [{^2,1 + "^2,3 ~ J2,t)(^1 + 

+ (■^2,4 - -^2,8 - 2 J2% + 2J2,15 - 2J2,ll)a2] . (6) 

Employing the calculated values of exchange parame- 
ters given (in K) in Fig. 4, together with calculated 
value of magnetic moment on iron /ii?e = 2.55/iB (cor- 
responding classical spin vector value 5=1.275), one 
can get from Eqs. (5) and (6) the Neel temperature 
T/v(ti)=217K, which is quite close to the experimen- 
tal value of T^^^ ~250K. Note that the account of ex- 
change bonds with adjacent ladders J only slightly in- 
fluences the value of Tjv (without the account of J we 
get Tjv(ri)=204K). However, the presence of J makes 
symmetry of exchange surrounding different for sublat- 
tices 1 and 2. This significantly influences the values 
and signs of exchange bonds along right and left sides 
of the spin ladder. 

Antiferromagnetic structure with spin configuration 
T2 is more complicated, however, using magnetic en- 
vironment symmetry one can establish the following 
equalities between average values of ^-projections of 
spins over all 16 sublattices (see Fig. 4, lower panel). 

C7i = — C75 = — CTio = 0-14, |J2 = —0-6 = — (Tg = (Tis, 
C73 = — C77 = —0-12 = 0-16, — (74 = (Tg = (Tn = — (Tl5-(7) 

Then following the same line of arguments as for calcu- 
lation of T/v(ti) (Eqs. (4-6)), one can get four equations 
for CTi , G2, (73 and (74 . Again using the calculated value 
of iron local magnetic moment jipe = "^ASfiB (5=1.215) 
and exchange integral values (in K) (see Fig. 4, lower 
panel), one can obtain Tjv(r2)=186K. This value is 36K 
lower than T/v(ti). This tells us that spin configuration 
Ti is, in fact, more favorable. It is probable that one of 



the reasons for this result is the presence of strong an- 
tiferromagnetic exchange 2 = ■^5,6="54K in T2 struc- 
ture. 

4. CONCLUSION 

In this work we have studied the band struc- 
ture, density of states and Fermi surfaces of recently 
discovered possible Fe-based superconductor BaFe2Se3 
(Bal23). In contrast to typical representatives of iron 
pnictide and chalcogenide systems [10, 11] and despite 
rather similar chemical composition, Bal23 has com- 
pletely different crystal structure. Instead of planes with 
square lattice of Fe-ions, Bal23 is a spin-ladder system 
with two-leg ladders of Fe ions, which form checker- 
board pattern in the plane perpendicular to the lad- 
der direction. Thus, not unexpectedly, the band disper- 
sions and Fermi surfaces are quite different from those in 
BaFe2As2 and (K,Cs)Fe2Se2. Complicated band struc- 
ture of Bal23 in the vicinity of the Fermi level illus- 
trates the importance of experiments on different kinds 
of doping in search for possible superconductivity in this 
system. 

For the undopcd Bal23 spin-ladder system, there is 
a problem of correct choice among two possible anti- 
ferromagnetic spin configurations, not easily discernible 
by neutron scattering ( "plaquettes" versus "zigzags"). 
We could distinguish these configurations by calcu- 
lating the Neel temperature for classical Heisenberg 
model, within the molecular field approximation. Corre- 
sponding Heisenberg exchange integral values were ob- 
tained from LSDA calculations. From our estimates 
we conclude that more favorable is "plaquette" config- 
uration with (calculated) T7v(ti)=217K (in contrast to 
"zigzags" with calculated rAr(T2)=186K). This conclu- 
sion agrees with experimental evidence of Ref. [7] and 
total LSDA energy calculations of Ref. [9] . 

This work is partly supported by RFBR grant 11- 
02-00147 and was performed within the framework of 
the Program of fundamental research of the Riissian 
Academy of Sciences (RAS) "Quantum mesoscopic and 
disordered structures" (12-0-2-1002). 
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